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A  Theoretical  Investigation  of  Laser-Sustained  Plasma  Thruster 
San-Mou  Jeng’  and  Dennis  Keefer** 

Center  for  Laser  Applications 
University  of  Tennessee  Space  Institute 
Tullahoma,  Tennessee  37388 
(615)  455-0631  ext.  338 


ABSTRACT 

A  numerical  code  has  been  successfully  devel¬ 
oped  for  the  investigation  of  thruster  performance 
using  a  laser-sustained  hydrogen  plasma  as  the  pro¬ 
pellant.  The  plasma  was  sustained  using  a  10.6/zm 
C02  laser  beam  which  is  focused  at  different  posi¬ 
tions  within  the  thruster.  The  physical  model  as¬ 
sumed  that  plasma  is  in  thermodynamically  equilib¬ 
rium  (LTE),  and  geometric  ray  tracing  was  adopted 
to  describe  the  laser  beam.  The  steady-state,  ax- 
isymmetric,  Navier-Stokes  equations  coupled  with  the 
laser  power  absorption  process  have  been  solved  nu¬ 
merically.  A  pressure  based  Navier-Stokes  numeri¬ 
cal  solver  using  body-fitted  coordinates  was  used  to 
calculate  the  laser-supported  rocket  flow  which  in¬ 
cludes  both  subsonic  and  supersonic  flow  regions. 
From  the  limited  parametric  study,  which  did  not 
try  to  optimize  the  rocket  performance,  it  was  found 
that  better  performance  was  obtained  when  the 
laser  beam  was  focused  closer  to  the  rocket  throat. 


NOMENCLATURE 

A  coefficient  in  difference  equation  for 

variable  <f> 

I  laser  ray  intensity 

J  Jacobian  of  inverse  coordinate 

transformation 
k  conductivity 

p  pressure 

91.92.  <73  geometric  relations  between  coordinate 
systems 

q  emission  coefficient  for  optically  thin 

radiation 
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R,S,T  source  term  in  equations  6,  7,  15, 
respectively 

a  distance  along  laser  ray 

U,  V  velocity  component  along  f  —  and 

T)— axis,  respectively 

u,  v  velocity  component  along  x—  and 

r— axis,  respectively 

f,  r?  sixes  of  curvilinear  coordinate  system 

<f>  general  dependent  variable 

p  density 

r  effective  diffusion  coefficient 

pi  viscosity 

a  absorption  coefficient 

Subscript 

E,  Wt  N,  S four  adjacent  nodes  to  P 

e,  w,  n,  s  four  surfaces  of  the  control  volume 

centered  at  P 
i  laser  ray  t 

P  nodal  point  to  be  solved  in  difference 

equation 
rad  radiation 

f ,  ri  partial  derivative  with  respect  to  f ,  rj 

<f>  for  dependent  variable  4> 

Superscript 

*  values  based  on  given  pressure  field 

'  corrected  values  according  to  corrected 

pressure  field 

INTRODUCTION 

This  paper  describes  the  numerical  calculations 
for  a  rocket  using  a  laser-sustained  plasma  as  the  pro¬ 
pellant  for  a  propulsion  device.  The  plasma  is  used 
to  absorb  and  convert  the  power  from  a  remote  laser 
beam  into  the  propellant  which  is  expanded  through 
a  nozzle  to  generate  the  required  thrust  for  a  space 
vehicle.  This  kind  of  propulsion  system  is  called  laser- 
thermal  propulsion.  Several  different  modes  of  oper¬ 
ating  a  laser-supported  rocket  have  been  proposed;  for 
example,  pulsed  laser-supported  detonation  wave;1 
cw  laser-sustained  plasma  (LSP).2-4  In  this  paper, 
the  discussion  is  limited  to  the  numerical  calculations 
of  a  rocket  flow  using  a  cw  laser-sustained  pure  hy¬ 
drogen  plasma  as  the  propellant. 

A  schematic  of  a  proposed  laser-supported  rocket 
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thruster  is  shown  in  Figure  1.  Laser  power  is  ab¬ 
sorbed  by  and  sustains  a  high  temperature  plasma 
core  (about  15,000-20,000  K),  which  mixes  with  the 
cold  buffer  flow  of  propellant  downstream,  heating  the 
flow  to  the  maximum  temperature  that  the  chamber 
wall  material  can  tolerate.  Then,  the  flow  is  accel¬ 
erated  through  a  conventional  nozzle  to  produce  the 
required  thrust.  Since  the  laser  power  absorption  by  a 
plasma  is  through  inverse  bremsstrahlung,  any  type  of 
propellant  can  be  used,  and  for  propulsion  purposes, 
the  lowest  molecular  weight  propellant,  hydrogen,  is 
the  best  candidate.  The  major  advantages  of  this  de¬ 
vice,  compared  to  chemical  propulsion,  are  that  the 
lowest  molecular  weight  propellant  can  be  used,  and 
that  the  maximum  temperature  of  the  working  fluid  is 
not  limited  to  the  flame  temperature  of  a  combustion 
process.  The  combination  of  high  temperature  and  a 
low  molecular  weight  propellant  results  in  a  very  high 
expected  specific  impulse.  An  estimated  specific  im¬ 
pulse  for  laser-sustained  plasma  propulsion  using  hy¬ 
drogen  propellant  is  in  the  range  of  1000-1500  sec.,6 
which  is  about  three  times  that  of  the  most  advanced 
chemical  propulsion  system. 


The  modeling  of  the  LSP  for  propulsion  began 
with  the  one-dimensional  work  of  Raizer.0  Several 
other  studies7,8  followed  this  analysis  but  used  more 
sophisticated  analytical  and  computational  methods. 
Unfortunately,  one-dimensional  models  have  an  inher¬ 
ent  weakness.  The  LSP  tends  to  be  highly  dependent 
on  laser  beam  geometry,  radial  convection  and  radial 
diffusion  transport  of  momentum  and  energy  in  the 
flow.  Keefer  and  coworkers9'10  developed  a  simplified 
semi-two-dimensional  model  which  considered  the  ra¬ 
dial  energy  diffusion  but  assumed  that  the  laser  beam 
and  the  streamlines  were  parallel.  Later,  dumb  and 
Krier11  extended  that  model  by  using  a  converging 
laser  beam,  as  well  as  using  more  realistic  gas  prop¬ 
erties.  However,  their  model  still  did  not  consider  the 
momentum  equations  but,  instead,  assumed  that  the 
radial  velocity  was  zero. 

Merkle  and  coworkers12,13  used  a  contemporary 
numerical  method  to  attack  the  full  two-dimensional 
problem.  They  succeeded  in  the  calculation  of  a  low 
temperature  LSP  (peak  temperature  about  4,000  K) 
in  a  low-mach  number  hydrogen  flow  seeded  with  ce¬ 
sium.  It  is  feasible  to  extend  their  approach  to  a  pure 
hydrogen  plasma,  but  the  higher  peak  temperature 
(about  16,000  K)  and  much  stronger  nonlinear  opti¬ 
cal  properties  of  the  pure  hydrogen  plasma  can  result 
in  numerical  stability  problems. 

Recently,  Jeng  and  Keefer14,16  obtained  the 
numerical  solutions  for  a  rigorous  two-dimensional 
model  with  realistic  gas  properties  for  the  LSP  in 
pure  argon  and  pure  hydrogen  flows.  Their  results  for 
the  argon  flows  agree  well  with  experimental  data10,17 
in  the  prediction  of  plasma  position,  size,  shape  and 


peak  temperature.  Jeng  et  al15, 18,19  have  also  stud¬ 
ied  the  effect  of  laser  power,  flow  rate,  static  pressure 
and  optical  geometry  on  the  LSP  flows.  Their  work 
also  indicated  that  there  are  no  fatal  errors  in  the 
laser  propulsion  concept,  and  the  critical  parameters 
of  plasma  size,  position,  radiation  heat  loss  and  en¬ 
ergy  conversion  efficiency,  can  be  controlled  using  an 
appropriate  optical  arrangement  and  a  suitable  flow 
configuration. 


The  objectives  of  the  current  study  are  to  de¬ 
velop  a  new  numerical  code  based  on  extended  plasma 
model  developed  by  Jeng  and  Keefer14,16  for  the  real 
nozzle  flow  which  includes  both  low  and  high  mach 
number  flows,  and  to  conduct  a  preliminary  paramet¬ 
ric  study  of  laser  rocket  performance.  There  was  no 
attempt,  in  this  paper,  to  determine  the  actual  rocket 
performance.  In  the  following  section,  the  theoretical 
method  used  in  the  study  is  described.  This  is  fol¬ 
lowed  by  a  discussion  of  model  calculations  and  some 
initial  parametric  studies. 


THEORETICAL  MODEL 


LASER  BEAM.  Geometric  optics  were  used  to 
describe  the  laser  beam  which  was  assumed  to  be  di¬ 
vided  into  a  finite  number  of  individual  rays.  Each 
ray  path  through  the  optical  system  to  the  plasma 
flow  was  calculated  using  real-ray  tracing.  Diffrac¬ 
tion  of  the  laser  beam  due  to  the  finite  aperture  of  the 
lens,  and  effects  due  to  the  inhomogeneous  refractive 
index  within  the  plasma  were  neglected  in  the  calcu¬ 
lation.  The  local  intensity  of  each  individual  ray  t  was 
described  by  Beer’s  law  written  in  the  the  following 
form. 


=  -a/. 


where  s,  is  the  distance  along  the  laser  propagation 
direction  for  ray  i,  and  a  is  absorption  coefficient  at 
laser  wavelength  and  depends  on  local  plasma  tem¬ 
perature  and  pressure. 


PLASMA  FLOW.  The  model  assumed  that  the 
flow  is  laminar,  steady-state,  axisymmetric  and  has 
variable  physical  properties.  The  pressure  of  the 
plasma  is  relatively  high  so  local  thermodynamic  equi¬ 
librium  (LTE)  for  the  plasma  can  be  assumed  with 
little  error.  Therefore,  the  plasma  can  be  described 
by  a  single  temperature,  and  its  physical  properties 
are  only  a  function  of  this  temperature  and  pres¬ 
sure.  The  model  also  considered  the  thermal  radia¬ 
tion  emitted  from  the  plasma  which  included  the  con¬ 
tributions  from  free  to  free  to  bound,  and  bound 
to  bound  transitions.  These  radiation  properties  are  a 
strong  function  of  wavelength,  and  the  current  model 
adopted  the  simplified  approach,  developed  by  Kemp 
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et  al20  and  Caledonia  et  al.21  to  solve  the  hydro¬ 
gen  plasma  radiation  transport  equations.  Radiation 
transport  was  divided  into  two  parts  depending  on 
wavelength:  for  wavelengths  longer  than  0.095  /im, 
the  optical  depth  is  long  and  an  optically  thin  ap¬ 
proximation,  in  which  all  radiation  escapes  from  the 
plasma,  was  used;  and  in  other  wavelength  regions, 
the  optical  depth  is  short  and  an  optically  thick  ap¬ 
proximation,  which  can  be  described  as  a  diffusion 
process,  was  assumed.  Following  the  above  assump¬ 
tions,  tae  equations  of  conservation  of  mass,  momen¬ 
tum  and  energy  for  the  flow  can  be  written  as  : 


d(P“)  1  djrpv) 

dx  r  dr 


(2) 


d(putt)  1  d(rpuv) 

dx  r  dr 


dp  d  /  9u\ 
dx  dx  \  dx  ) 

1  d  (  9u\  d  (  9u\  1  d  (  dv\ 

+ rA'ai)  *-,^1;) 
_ £.  [20  .  i  a(™)M 

dx  [  3  \9x  r  dr  )  J  ’ 


(3) 


where  kejf  is  the  sum  of  molecular  thermal  conduc¬ 
tivity  and  radiation-induced  thermal  conductivity  in 
the  optically  thick  limit;  qrad  is  radiative  heat  loss 
in  the  optically  thin  limit;  aE/,-  represents  the  local 
power  addition  from  the  laser  beam  to  the  flow;  and 
as  mentioned  before,  the  physical  properties  involved 
in  equations  (2-5)  are  a  function  of  the  local  temper¬ 
atures  and  pressures.  Reference  15  has  a  detailed  de¬ 
scription  of  the  thermophysical  and  optical  properties 
sources  used  in  this  work. 

In  the  calculations,  adiabatic,  nonslip  and  zero 
emissivity  conditions  were  assigned  at  the  thruster 
wall.  Total  pressure,  total  temperature  and  zero  ra¬ 
dial  velocity  were  specified  at  the  inlet,  and  upstream 
running  Riemann  invariants  extrapolation  plus  the 
isentropic  relations  were  used  to  calculate  other  prop¬ 
erties.  At  the  exit,  the  properties  were  linearly  ex¬ 
trapolated  from  the  interior  flow  domain. 

METHOD  OF  COMPUTATION 

COORDINATE  TRANSFORMATION.  The  set  ■ 
of  transport  equations  (3-5)  can  be  written  in  the 
cylindrical  system  coordinates  for  the  dependent  vari¬ 
able,  <f>,  in  the  following  generalized  form: 


d(puv)  1  d(rpvv) 
dx  +  r  dr 


dp  d  ( 

--Tr  +  TxV 
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rdr 
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( ( dv\  /<9u\  \  A  dudv  dv  du 
+HV^)  +{d~r)  )  ^VYxTr^^YxTr 

1  / 4  v2  4  f  dv  dv\ \ 
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5^(W)  +  A( fW)- A(r*rg) 

9  .  a,  d<f>. 

+r.R*{x,r),  (6) 

where  T*efj  is  the  effective  diffusion  coefficient  and 
R*(x,  r)  is  the  source  term.  The  source  term  for  the 
momentum  equations  (3,4)  also  includes  the  diffu¬ 
sion  terms  from  the  dilation  effect.  When  new  in¬ 
dependent  variables  f  and  rj  are  introduced,  Eq.  (6) 
changes  according  to  the  general  transformation  f  and 
rj.  A  schematic  illustration  of  the  relations  between 
the  physical  domain  and  the  transformed  domain  is 
shown  in  Figure  2.  Equation  (6)  is  rewritten  in  f  and 
T]  coordinates  as  follows: 


19,  ....  13,  , 

7^(r^)+7-(rp^)  = 


1  9_ 

Jd£ 


Cr  *r 


r  9*7 

J  +  J  dn 


+r-S*U,r,) 


r*r 

(1) 


U  =  urn  -  vxn, 


where 


(3) 


qi  =  x2n+r2,  (10) 

q2  =  x6Xr,  +  r(rn,  (11) 

qz  =  x]+y2i,  (12) 

J  =  ~  xnri}  (13) 

and  S*(f,  rj)  is  the  source  term  in  the  f  and  ij  coor¬ 
dinates. 

NUMERICAL  .  GORITHM  A  nonstaggered 
grid  system,  as  shown  in  Figure  2,  is  adopted  in  the 
calculations.  For  a  typical  grid  node  P  are  enclosed 
in  its  cell  and  surrounded  by  its  neighbors  N,  S,  E, 
and  W.  The  integral  of  Eq.  7  over  its  control  volume, 
by  arbitrary  taking  Af  =  1  and  A  r]  =  1,  becomes 


92  =  ZfXrj  +  r(rn, 
93  =  x\  +  y|. 


(rpU 4>)e  -  ( rpU<t>)w  +  ( rpV4>)n  -  (rpV<f>), 

rr*  rrr* 

=  ~j~[qi<r>(  —  q^t))  -  —j-{qi<t>i-q24>^) 

■  "  e  J  w 

rT*  rT*' 

+  —j-(-q2<t>i  +  93$rj)  -  -J-(  — 92^«  +  93^f,) 

-  n  J  $ 

+r  ■  J  ■  (14) 

Second-order  center  differencing  is  used  to  dis¬ 
cretize  the  above  equation,  except  for  the  convective 
terms  when  the  cell  Reynolds  number  is  greater  than 
2.  In  that  case,  first  order  upwind  differencing  is  used 
for  the  convective  terms.  A  staggered  grid  system,  for 
body-fitted  coordinates,  has  been  widely  used,22  and 
it  has  been  known  that  an  oscillatory  field  will  be  pro¬ 
duced  using  the  nonstaggered  grid  system.  The  arti¬ 
ficial  viscosity  term  proposed  by  Rhie  and  Chow23 
to  cure  the  problem  was  also  implemented  in  this 
study.  The  cross  derivatives  in  the  diffusion  terms 
of  Eq.  (14)  are  the  results  of  the  nonorthogonal  coor¬ 
dinate  system  and  are  usually  small.  In  order  to  use 
the  solution  procedure  for  the  five-point  unknowns, 
these  terms  are  combined  with  the  source  term  and 
treated  as  known  quantities  and  the  resulting  rela¬ 
tion  between  4>p  and  the  neighboring  values  can  be 
written  as  : 


Ap<j>p  =  Ae4>e  +  Aw4>w  +  An4>n  +  As<f>s 

(15) 

where  the  coefficients  A  involve  the  flow  properties  of 
convection,  diffusion,  area,  etc.,  and  rj)  includes 
the  source  term  rj)  and  the  cross  derivatives  in 

the  diffusion  terms. 

PRESSURE  CORRECTION  PROCEDURE. 
While  the  general  transport  equation,  Eq.  (15),  was 


used  to  solved  momentum  and  energy  equations,  the 
pressure  field  must  be  solved  using  the  continuity 
equation.  The  procedure  used  in  this  work  is  an  adop¬ 
tion  of  algorithm  developed  by  Rhie,24  and  only  a 
brief  discussion  is  presented  in  this  section.  The  pres¬ 
sure  equation  used  in  the  adopted  procedure  is  similar 
to  that  of  SIMPLER25  and  PISO,20  27  except  that  the 
density  term  is  treated  implicitly  for  the  compressible 
flow.  The  basic  idea  is  to  use  the  following  relations 
to  correct  the  mass  flux  obtained  from  momentum 
equations 


pu  =  p*u*  +  p'u*  +  p*u'  (16) 

pv  =  p*t>*  +  p'u*  +  p*t/  (17) 

and  the  discretized  pressure  equation  becomes  a 
transport  equation  for  compressible  flow  instead  of 
Poisson  equation.  The  ratio  of  diffusion  term  to  con¬ 
vection  term  is  dependent  on  Mach  number,  and  as 
Mach  number  increases,  the  convection  term  becomes 
dominant.  The  split  operator  concept  similar  to  PISO 
was  also  adopted  in  the  current  study  to  solve  the 
pressure  correction  equation. 

RESULTS  AND  DISUCSSION 


Figure  3  shows  the  thruster  dimensions  used  in 
the  calculations.  It  includes  a  cylindrical  section  up¬ 
stream  of  the  convergent-divergent  nozzle.  The  throat 
and  cylindrical  diameters  are  6  and  40.74  mm  respec¬ 
tively,  and  the  exit  to  throat  area  ratio  is  6.  The  coor¬ 
dinate  origin  used  to  present  the  calculated  results  is 
at  the  center  of  the  throat.  Although  Figure  1  shows 
that  the  hydrogen  flows  into  the  chamber  radially, 
in  the  calculations  reported  here  the  inlet  flow  direc¬ 
tion  was  assumed  axial.  The  total  pressure  and  total 
temperature  were  specified  at  the  inlet.  In  the  tradi¬ 
tional  chemical  rocket,  the  total  pressure  has  a  minor 
effect  on  the  expected  value  of  the  specific  impulse, 
but  in  the  laser-supported  rocket,  the  pressure  has  a 
substantial  influence  on  the  absorption  and  radiation 
processes  that  directly  affect  the  available  enthalpy  of 
the  propellant  and  the  specific  impulse  of  the  thruster. 
Based  on  our  previous  parametric  studies15  for  laser- 
sustained  hydrogen  plasmas  in  a  constant  area  pipe 
flow,  it  was  found  that  3  atm  total  pressure  was  ad¬ 
equate  for  the  thruster  upstream  condition,  and  this 
value  was  used  in  the  calculations.  The  total  tem¬ 
perature  at  the  inlet  was  specified  as  1500K  in  the 
calculations,  based  on  the  expected  use  of  a  regenera¬ 
tive  cycle  to  recover  the  optically  thin  radiation  from 
plasma. 

The  focusing  lens  used  in  the  calculations  has  a 
nominal  focal  length  203  mm,  and  References  14  and 
15  has  detailed  discussion  of  this  lens  geometry.  The 
actual  dimensions  of  this  lens  were  adopted  in  the  ray 
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tracing  for  the  laser  beam,  so  the  spherical  aberration 
of  the  lens  was  also  included  in  the  calculations.  The 
laser  was  assumed  to  have  10.6  /zm  wavelength  and  a 
collimated,  40  mm  diameter  Gaussian  beam  incident 
on  the  lens.  The  only  preliminary  parametric  study 
of  the  rocket  performance  reported  in  this  paper  was 
to  focus  the  laser  beam  at  different  axial  locations 
using  the  same  lens  arrangement.  Table  1  lists  the 
operating  conditions,  including  one  case  without  laser 
power  addition. 

In  the  calculations,  nonuniformly  spaced  grids  in 
the  physical  space  were  transformed  into  uniformly 
spaced  grids  in  both  radial  and  axial  directions  in  the 
transformed  space,  as  illustrated  in  Figure  4.  From 
previous  studies  14-15-18-19  cf  LSP  within  a  constant 
diameter  pipe,  the  position  of  the  high  temperature 
plasma  zone  was  found  to  be  highly  dependent  on  the 
flow  configurations  and  the  optical  arrangements,  and 
it  was  not  easy  to  prelocate  the  region  in  the  physical 
plane  where  a  dense  grid  distribution  was  needed  to 
accommodate  the  high  temperature  gradients.  There¬ 
fore,  a  preliminary  calculation  was  first  performed  for 
each  case  to  locate  the  high  temperature  zone,  then 
the  grid  distribution  was  refined  and  used  for  the  final 
calculation. 

The  newly  developed  computer  code  was  first 
used  to  calculate  the  thruster  flow  without  consid¬ 
ering  laser  heat  addition,  and  the  results  were  used 
to  check  the  behavior  of  the  computer  program  and 
to  serve  as  benchmark  to  compare  with  the  other 
calculated  cases  which  included  the  laser  heat  ad¬ 
dition.  The  Mach  number  contours  for  case  no.  1 
are  shown  in  Figure  5.  Because  the  nozzle  expan¬ 
sion  angle  is  only  6°,  the  contours  outside  of  the  wall 
where  diffusion  transport  is  not  important,  are  almost 
straight,  and  are  similar  to  one-dimensional  nozzle 
calculations.  Since  the  nonslip  condition  was  applied 
in  the  calculations,  the  Mach  number  near  the  wall  is 
smaller  than  that  near  the  centerline  at  the  same  axial 
location.  This  demonstrates  that  the  newly  developed 
computer  program  can  also  handle  the  boundary  layer 
in  supersonic  flow.  In  this  test  case,  there  were  only  a 
few  grids  within  the  boundary,  and  we  made  no  spe¬ 
cial  attempt  to  resolve  the  boundary  layer  to  good 
numerical  accuracy.  If  needed,  more  grids  can  be  in¬ 
cluded  in  that  region  and  a  more  accurate  solution 
can  be  obtained. 

Figure  6  shows  isotherms  for  case  no.  2,  and  the 
outer  laser  beam  boundary  was  also  plotted  in  the 
figure  as  a  dashed  line.  The  plasma  core  is  located 
upstream  of  the  laser  focus,  and  the  distance  between 
the  laser  focus  and  the  high  temperature  plasma  core, 
which  depends  on  flow  configurations  and  optical  ar¬ 
rangements  as  discussed  in  the  early  section  of  this 
paper,  is  10  mm.  Downstream  of  the  laser  focus,  the 
mixing  process  quickly  reduces  the  centerline  temper¬ 


ature.  In  the  throat  region,  only  a  small  portion  of 
the  hydrogen  flow  near  the  axis  is  heated,  and  a  major 
portion  of  the  flow  is  still  unheated.  For  an  optimized 
rocket,  one  of  the  design  criteria  is  to  have  a  uniformly 
distributed  temperature  at  the  throat  in  order  to  re¬ 
duce  the  two-dimensional  losses.  In  the  calculated 
case,  this  condition  was  not  satisfied.  Near  the  wall 
of  the  nozzle  exit,  there  is  a  local  hot  spot  which  b  due 
to  the  viscous  dissipation  effect  within  the  boundary 
layer. 

Laser-sustained  plasmas,  operated  at  30  kW  laser 
power,  absorb  almost  all  of  the  incoming  laser  beam 
power.  In  the  calculated  case  no.  2,  only  397  W 
of  laser  power  was  transmitted  through  the  plasma 
which  either  reached  the  thruster  wall  or  was  trans¬ 
mitted  through  the  throat.  The  wall  material  could 
easily  withstand  this  small  level  of  transmitted  power 
density.  The  optically  thin  radiation  from  the  plasma 
is  also  an  important  factor  for  a  laser-supported 
rocket  design.  It  represents  the  portion  of  laser  power 
that  is  not  directly  converted  into  the  flow  power. 

However,  this  radiation  power,  like  a  part  of  the  trans¬ 
mitted  laser  power,  will  be  absorbed  by  the  wall,  and 
by  using  a  regenerative  cycle,  part  of  this  energy  can 
be  recovered  and  used  to  heat  the  hydrogen  from  its 
storage  temperature  up  to  the  propellant  inlet  tem¬ 
perature  of  the  thruster.  For  case  no.  2,  the  optically  j 

thin  radiation  was  17.9  kW  which  was  only  30  %  of  i 

the  energy  needed  to  heat  hydrogen  to  the  assigned  j 

inlet  temperature,  1500  K. 

Pressure  and  Mach  number  contours  for  case  no. 

2  axe  shown  in  Figure.  7  and  Figure.  8,  respectively. 

Pressure  contours  downstream  of  the  laser  focus  in  1 

Figure  7  are  similar  to  the  calculated  results  for  case 
no.  1.  Those  have  only  a  very  small  radial  variation  I 

and  imply  that  the  nonuniformly  distributed  temper¬ 
ature  field  has  little  effect  on  the  pressure  field.  How-  I 

ever,  a  high  pressure  zone,  even  higher  than  the  inlet 
static  pressure  of  the  thruster  (3xlOsIV/m2),  occurs 
ahead  the  high  temperature  plasma  zone.  The  Mach 
number  at  the  inlet  is  about  0.019  and  the  maximum 
calculated  Mach  number  is  2.605  at  the  nozzle  exit. 

The  Mach  number  contours  in  the  supersonic  region 
are  very  complicated,  compared  to  that  of  case.  no. 

1  in  Figure  5,  although  the  velocity  monotonically  in¬ 
creased  from  the  wall  to  the  centerline.  The  Mach 
number  at  the  centerline  is  always  smaller  than  away 
from  the  centerline,  since  the  temperature,  as  well  as 
sonic  speed,  is  higher  at  the  centerline.  The  unheated 
region  still  has  very  straight  Mach  number  contours 
similar  to  those  shown  in  Figure  5  for  case  no.  1.  The 
mixing  region  between  hot  center  flow  and  the  colder 
unheated  flow  shows  very  complicated  contours.  Mix¬ 
ing  in  that  region,  which  implies  variations  of  sonic 
speed  due  to  different  temperature,  molecular  weight, 
and  specific  heat,  are  the  major  factors.  The  Mach 
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number  near  the  wall  is  smallest  due  to  the  bound¬ 
ary  effect,  and  the  Mach  number  on  the  centerline  is 
always  smaller  than  that  away  from  the  centerline  in 
the  subsonic  portion  of  the  flow. 

In  the  previous  studies15,16,  it  was  found  that 
laser  heat  addition  to  the  flow  will  generate  a  very 
complicated  velocity  field,  even  in  a  constant  area 
pipe.  The  same  behavior  was  also  observed  in  the 
current  results.  Figure  9  shows  the  radial  velocity 
contours  for  case  no.  2.  Near  the  high  temperature 
plasma  core  region,  both  positive  and  negative  radial 
velocity  occurs.  In  the  supersonic  region,  a  local  min¬ 
imum  exists  in  the  mixing  zone  between  heated  and 
unheated  flow  for  case  no.  2,  and  the  radial  veloc¬ 
ity  increases  from  the  centerline  up  to  the  edge  of 
the  boundary  for  the  unheated  case  no.  1.  Nonuni- 
formly  distributed  temperature  at  the  throat  is  the 
main  cause  for  the  above  results. 

Figure  10(a)  and  10(b)  show  isotherms  for  cases 
no. 4  and  no.  5,  respectively.  Their  laser  focuses  are 
closer  to  the  throat  than  case  no.  2,  and  the  focus  is 
within  the  nozzle  contracting  section  for  case  no.  5. 
As  the  laser  focus  moves  closer  to  the  throat,  the  size 
of  the  plasma  becomes  smaller.  Since  the  optically 
thin  radiation  is  proportional  to  the  plasma  size,  the 
radiation  loss  for  case  no.  5,  9.9  kW,  is  much  smaller 
than  15.9  kW  and  18.6  kW  for  cases  no.4  and  no. 
2,  respectively.  Because  the  transmitted  laser  power 
is  small  and  similar  for  all  cases,  the  energy  directly 
converted  to  the  flow  is  highest  for  case  no.  5.  The 
laser  beam  also  heats  a  larger  portion  of  the  flow  at 
the  throat  in  case  no.  5  than  other  cases,  and  this 
produces  better  thruster  performance. 

Table  2  lists  the  calculated  properties  for  the 
five  test  cases.  The  optically  thin  radiation  heat  loss 
ranges  from  61  %  to  33  %  of  the  laser  power,  and  the 
transmitted  laser  power  is  relatively  small  by  com¬ 
parison.  The  peak  axial  velocity,  which  is  always  lo¬ 
cated  on  the  centerline  at  the  exit,  increases  from  4763 
m/sec  for  the  unheated  thruster  to  18957  m/sec  for 
case  no.  5.  The  maximum  temperature  of  the  plasma 
shows  little  dependence  on  the  variation  of  laser  focal 
position  among  the  test  cases.  The  peak  Mach  num¬ 
ber  at  the  exit  plane  is  located  within  the  unheated 
flow  for  all  test  cases,  and  it  increases  as  the  laser 
focal  point  moves  downstream.  The  mass  flow  rate 
decreases  as  the  laser  focal  point  moves  toward  the 
throat. 

SUMMARY  AND  CONCLUSIONS 

The  flow  within  a  laser-sustained  plasma  rocket 
thruster  has  been  successfully  calculated.  The  physi¬ 
cal  model  was  based  on  steady-state  two-dimensional 
Navier-Stokes  equations  for  the  flow  and  geometric 
real-ray  tracing  for  the  laser  beam.  Realistic  thermo¬ 


physic  and  optical  properties  based  on  locally  ther¬ 
modynamical  equilibrium  assumptions  were  incorpo¬ 
rated  into  the  calculations,  and  radiation  heat  trans¬ 
fer  from  the  plasma  was  simplified  using  optically  thin 
and  optically  thick  assumptions.  The  governing  equa¬ 
tions  were  first  transformed  into  body-fitted  coordi¬ 
nates  and  were  solved  numerically.  A  pressure-based 
Navier-Stokes  solver  with  a  split  operator  concept  was 
used  for  the  numerical  procedure. 

Several  laser  propulsion  cases  were  calculated  us¬ 
ing  the  newly  developed  computer  code.  The  noz¬ 
zle  geometry,  flow  configuration  and  optical  arrange¬ 
ments  were  not  optimized  for  rocket  performance. 
However,  some  useful  information  was  obtained  from 
these  preliminary  test  cases.  The  performance  of  the 
rocket  is  highly  dependent  on  the  laser  focal  position. 
The  maximum  velocity  is  larger  when  laser  focal  point 
moves  closer  to  the  throat.  Generally  speaking,  the 
laser  beam  only  heated  a  small  portion  of  the  hydro¬ 
gen  propellant,  and  the  best  thruster  performance  will 
not  be  achieved  due  to  the  two-dimensional  losses.  In 
the  future,  an  extensive  parametric  study  will  be  re¬ 
quired  to  optimize  the  rocket  performance,  and  one 
of  the  major  goals  of  the  study  is  to  design  the  rocket 
to  produce  a  more  uniformly  distributed  temperature 
at  the  throat. 
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Figure  1.  Schematic  of  a  laser-sustained 
plasma  propulsion  thruster  chamber. 
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Figure  2.  Finite-difference  grid  representation: 
(a)  physical  plane,  (b)  transformed  plane. 
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Figure  4.  Schematic  of  grid  transformation  (a) 
physical  plane,  (b)  transformed  plane. 
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Figure  5.  Mach  number  contours  for  case  num¬ 
ber  1. 
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Figure  8.  Mach  number  contours  for  case  num 
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Figure  6.  Temperature  contours  for  case  num¬ 
ber  2. 
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Figure  9.  Radial  velocity  contours  for  case 
number  2. 
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Figure  7.  Pressure  contours  for  case  number 
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‘  Table  1:  Test  Conditions 


Case  no. 

Laser  power  fkW)  focus  fmml 

1 

0 

2 

30 

-85 

3 

30 

-50 

4 

30 

-35 

5 

30 

-20 

(b) 


Figure  10.  Temperature  contours:  (a)  case 
number  4,  (b)  case  number  5. 


Table  2:  Calculated  properties  for  five  test  cases 


Case  no.  1 

2 

3 

4 

5 

Optically  thin.  rad.  (kW)  0. 

18.6 

17.9 

10.0 

9.9 

Transmitted  power  (W)  0. 

397 

395 

330 

320 

Max.  axial  vcl.  (m/sec)  4763 

12513 

13403 

15043 

18957 

Max.  temperature  (K)  1500 

16854 

10872 

10878 

17293 

Max.  mach  number  2.50 

2.u  i 

2.64 

2.68 

2.78 

Mass  flow  rate  fg/sec)  2.88 

O  QO 

2.74 

2.69 

2.00 
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